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ABSTRACT 

We investigate the extended source size effects on gravitational lensing in 
which a lens consists of a smooth potential and small mass clumps ("substruc- 
ture lensing"). We first consider a lens model that consists of a clump modeled 
as a singular isothermal sphere (SIS) and a primary lens modeled as an exter- 
nal background shear and convergence. For this simple model, we derive analytic 
formulae for (de)magnification of circularly symmetric top-hat sources with three 
types of parity for their lensed images, namely, positive, negative, and doubly 
negative parities. Provided that the source size is sufficiently larger than the 
Einstein radius of the SIS, we find that in the positive (doubly negative) par- 
ity case, an extended source is always magnified (demagnified) in comparison 
with the unperturbed macrolens system, whereas in the negative parity case, the 
(de) magnification effect, which depends on the sign of convergence minus unity 
is weaker than those in other parities. It is shown that a measurement of the 
distortion pattern in a multiply lensed image enables us to break the degeneracy 
between the lensing effects of clump mass and those of clump distance if lensing 
parameters of the relevant macrolens model are determined from the position 
and flux of multiple images. We also show that an actual density profile of a 
clump can be directly measured by analyzing the "fine structure" in a multiply 
lensed image within the Einstein radius of the clump. 

Subject headings: cosmology: theory - dark matter - gravitational lensing - 
large-scale structure of universe 

1. Introduction 

The standard cold dark matter (CDM) scenario predicts the presence of several hundred 
small-mass clumps or "subhalos" (< 1O 8 M ) in a galaxy-sized halo (~ 1O 12 M ), while the 
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observed number of dwarf galaxies around the Milky Way is only a dozen (Klypin et al. 
1999; Moore et al. 1999). This suggests the presence of several hundred subhalos holding 
few or no stars around such a galaxy. Gravitationally lensed QSOs with quadruple images 
have recently been used for putting a limit on the surface density and the mass of such 
invisible substructures (Mao & Schneider 1998; Metcalf & Madau 2001; Chiba 2002; Metcalf 
& Zhao 2002; Dalai & Kochanek 2002; Bradac et al. 2002). Statistical analyses on these 
quadruple QSO-galaxy lensing systems show that the (de) magnification of lensed images and 
its dependence on their parities can be explained by the presence of substructures along the 
line of sight to the images, in contrast to the lens models relying on "smooth perturbation" 
in the gravitational potential of a macrolens (Metcalf & Zhao 2002; Keeton et al. 2002; 
Evans & Witt 2003). 

However, there are two ambiguities in these QSO-galaxy lensing systems. First, the 
subhalo masses are poorly measured. We can put meaningful constraints on these masses by 
observing substructure lensing of extended sources with various sizes (Moustakas & Metcalf 
2003; Dobler & Keeton 2005). If the angular source size is sufficiently larger than the Einstein 
radius of the subhalo, then the images lensed by a smooth galaxy halo will not be perturbed 
by a subhalo at all. Thus, by using a source with a large angular size, we can put a stringent 
constraint on the lower limit of the subhalo mass. Various examples of an extended source 
are available in QSO-galaxy lensing systems. For instance, a QSO core component has a 
typical size of < 1 x (z//5GHz) _1 pc in the radio band (Kameno et al. 2000; Kadler et al., 
2003), which is to be compared with the Einstein radius of ~ 1 x (M/ (3 x 10 5 M Q )) 1/2 pc for 
a point-mass clump with mass M. In addition, a surrounding hot dust torus around a QSO 
nucleus is supposed to be as large as 0.1 — 1 pc, and the blackbody emission from this hot 
dust is observable in the infrared band (e.g., Agol, Jones, & Blaes 2000; Chiba et al. 2005). 
Cold dust components at temperatures of 20 — 50K around a QSO nucleus can be as large 
as several kpc (e.g., Puget et al. 1996). 

Second, the locations of subhalos along the line of sight are poorly constrained. In 
fact, extragalactic substructures other than those associated with primary lensing halos 
can also contribute significantly to QSO-galaxy lens systems. Using iV-body simulations, 
Wambsganss et al. (2004) showed that the fraction of cases for which more than one lens 
plane contributes significantly to the multiply imaged system is very large for a source at 
high redshift z > 3. Metcalf (2005) showed that all the cusp caustic lens anomalies can 
also be explained by extragalactic ACDM halos with a mass range of 10 s — 1O 9 M . In 
order to determine the distances to these substructures, we need to have more information 
in addition to the position and flux of multiple images. Future high-resolution mapping of 
extended images in a QSO-galaxy lensing system may provide us useful information about 
the masses, distances, and density profiles of such substructures (Inoue & Chiba 2003, 2005). 
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Such an observation can be achieved by next-generation space VLBI such as the VLBI Space 
Observatory Program 2 (VSOP2; Hirabayashi et al. 2001) and submillimeter interferometers 
such as the Atacama Large Submillimeter Array (ALMA). 

To address these issues, we need to understand substructure lensing effects for an ex- 
tended source. However, to date little attention has been paid to the effects of a source 
size that is comparable to or larger than the Einstein radius of a perturber although lensing 
of such an extended source by a point mass (i.e., without a background, primary lens) has 
been studied extensively (e.g., Witt & Mao 1994). In comparison with lensing systems in 
which the primary lens is a cluster of galaxies, a galaxy-sized halo can be represented by a 
relatively small number of parameters since it has many symmetries. This allows us to treat 
the lensing effect of small-mass substructures as a local perturbation to the macrolensing 
effect in QSO-galaxy lensing systems. For point sources, the dependence on image parities 
of their (de)magnification has been investigated (Finch et al. 2002; Keeton 2003). It is 
of great importance to know to what extent the finite source size affects such systematic 
(de) magnification in substructure lensing. 

In this paper, we investigate the extended source effects in substructure lensing. In 
section 2 and 3, we show that for an SIS lens (with or without background shear and 
convergence), the flux of a circular source can be significantly altered even if the source size 
is much larger than that of the Einstein ring. In section 2, we consider a single SIS model 
without background for a circularly symmetric top-hat or Gaussian surface brightness profile 
of a source. In section 3, we explore SIS lens systems with an external background shear 
7 and convergence k relevant to substructure lensing. We derive analytic (de) magnification 
formulae based on astrometric shifts for a circular symmetric top-hat source for three types of 
parity of an image, namely, positive, negative and doubly negative parities. We show that in 
the large source size limit, the (de) magnification effect depends systematically on the parity 
of an the image, and it is prominent even for sources somewhat larger than the Einstein ring. 
We also discuss the mean (de) magnificat ion perturbation for a point source. In section 4, 
we study a simple method of breaking the degeneracy between the lensing effects of subhalo 
mass and distance using astrometric shifts of a macrolensed image. In section 5, we explore 
a method to directly measure the mass density profile of each lens perturber. In section 6, 
we summarize our results. 



2. SIS lens without background 

In this section, we show that for an SIS, the magnification effect is still prominent even 
if the source size is larger than the Einstein radius. This implies that if the source size 
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dependence of magnification is measured, then one can determine the Einstein radius of an 
SIS lens if the source center is fixed. Note that an SIS is widely used in the literature for 
representing the density profile of CDM subhalos with 10 7 — 10 9 M Q relevant to substruc- 
ture lensing (Metcalf & Madau 2001; Dalai & Kochanek 2002), although it is somewhat 
more concentrated than the so-called Navarro-Frenk- White (NFW) profile (Navarro, et al. 
1996) near the central cusp. For simplicity, we ignore the effect of a background shear and 
convergence as well as the ellipticity of the lens, in this section. 

Let X = {0 x i,d X 2) and 6 y = (6 y i,0 y 2) be two-component vectors of angles in the sky 
representing the image position and the source position, respectively. The lens equation for 
an SIS is then 

Q 

y = X - 0E77T7' (!) 

where the Einstein radius 6e is written in terms of velocity dispersion a and light velocity c 

as 

where D LS and D s are angular diameter distances between the lens and the source, and 
between the observer and the source, respectively. In what follows, we set 6e = 1 without 
loss of generality. In terms of complex variables z = Q x \ + i9 x2 and ( = Q y \ + i0 y 2, the 
normalized lens equation is expressed as 

\z\ 

which has a set of solutions 

z+ = (l+ IC|)4 for V C, z. = (|C| - 1) ||p for |C| < 1. (4) 

corresponding to an image position with positive parity (denoted by the subscript plus sign) 
and that with negative parity (denoted by the subscript minus sign). 



2.1. Top-hat source 

In this subsection, we explore a circularly symmetric top-hat source with a radius L in 
units of the Einstein radius with a constant surface brightness, located at a distance £ from 
the center of an SIS lens. 

In terms of a complete elliptic integral of the first kind K(k) and that of the second 
kind E(k), the magnification factors of a circular top-hat source corresponding to either a 
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positive parity image z + or a negative parity image Z- can be written as 



^(iL + Co)E(k) + (L-Co)K(k)j , k = 




1/2 



li±(L,(o) = 1± 



(5) 



To understand the effect of the finite source radius L, first, we study the case in which 
a circular top-hat source is centered at the lens center (Co = 0). From equation (5), one 
can find that the magnification factor fi = fi_ diverges as /i ~ 4/L in the small size limit 
L — > 0, similar to that of a point mass lens with a unit Einstein radius, /i ~ 2/L (Witt & 
Mao 1994). However, in the large-size limit L — > oo, the magnification factor for an SIS 
lens converges more slowly to unity, as fi ~ ~ 1 + 2/L + l/L 2 , than that of a point 
mass lens, fx ~ 1 + 2/L 2 (Witt & Mao 1994). Therefore, for an SIS, the magnification effect 
cannot be negligible even if the source size is greater than the Einstein radius. For instance, 
the source flux can be altered by 20 percent even if the source size is 10 times larger than 
the Einstein radius. This is because the deviation of the photon path coming from regions 
outside the Einstein radius is larger for an SIS lens than that for a point mass lens with the 
same Einstein radius. This behavior is reasonable, since the deflection angle for an SIS lens 
is constant with increasing \0 X \, whereas it decays as 1/6 X for a point mass lens. 

Next, we study the case in which a circular top-hat source is placed at an off-center 
position (0, Co 7^ 0) in the lens plane. For L > land(o > 1, the magnification factor can be 
approximately given by (see Appendix A) 



where s > controls smoothness of 5fi around L = Co- As shown in Figure 1, the magni- 
fication factor ii keeps its value fi ~ 1 + 1/Co until the boundary of the source touches the 
edge of an Einstein ring. Then it starts to increase until the source completely includes the 
Einstein ring adding an extra contribution 5fi = l/L 2 corresponding to the "second" image, 
which is shown in the top right and bottom left panels in Figure 2. As the source radius L 
increases, fx gradually decreases as ii ~ ii + + l/L 2 . Because /i + ~ 1 + 2/L for L 3> 1, the 
magnification factor can again be significantly altered even if the source size is larger than 
that of the Einstein ring provided that the Einstein ring is totally "occulted" by the source. 



K L > Co) ~ V+(L, Co) + 5fi(s, L, Co), L > 1, Co > 1. 



(6) 



where 



Ms,£,Co) 



Erf(s(L - Co)) ~ Erf(-s(L - Co)) + 2 
AL 2 



(7) 
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2.2. Gaussian source 



In this subsection, we explore a circularly symmetric Gaussian source located at distance 
Co from the center of an SIS lens. The surface brightness of a Gaussian circular source with 
standard deviation L is 

icr 



2nL 2 



cxp 



2L 2 



The magnification factor for a Gaussian circular source is then 



all 



(8) 



(9) 



For an image with positive parity, equation (9) can be reduced to 



^ r-co r-2-K dim Z 

—— / dR / d(f) Re z + — — — -Rexp 
2vrL J J 



d(t> 



R 2 



2L 2 



(10) 



Expanding equation (10) in Co, we obtain the magnification factor /x in the large source size 
limit, L — > oo, as 



G 



1 + (tt/2) 1 / 2 L- 1 + L- 2 /2. 



(11) 



As one can see in Figure 3, the contribution from an image with negative parity is negligible 
in the large source size limit. As in the top-hat model, the order of the leading term in /i — 1 
is 0(L~ V ). Therefore, the magnification effect for a Gaussian source cannot be negligible 
even if the source size is greater than the Einstein ring. 

To estimate \i in the small source size limit L <C 1, we introduce a cutoff radius R cut ~ L 
such that the contribution of integrand in equation (10) for R > R cnt is negligible. If we 
further assume that R cut <C 1 < Co, then equation (10) yields 



,,G -"-cut 



-5 



(12) 



In the limit Co ~~ * °°! 110 magnification is expected, i.e., fj, — 1. Therefore, we obtain the 
same formula for the small source size limit, /i G ~ 1 + Co" 1 f° r L <^ 1 and Co > 1, as in the 
top-hat model. 

As one can see in Figure 3, for L ^> Co + 1 the behavior of the magnification factor for 
a Gaussian source [P as a function of L is similar to that of \i for a top-hat source, if an 
appropriate scaling for the source radius is carried out. Therefore, for an SIS, we conclude 
that the surface brightness profile does not play an important role in magnification in the 
large source size limit, provided that the source has a circular symmetry. 
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However, in the Gaussian model, distorted images have two distinct features depending 
on the position of the SIS : (1) a pair of bright and dark spots (dipole structure) at the 
position of the SIS where the spatial gradient in the surface brightness V/ does not vanish 
(Figure 4, top), and (2) a dark spot (monopole structure) at the position of the SIS at the 
peak in the surface brightness where V/ = (Figure 4, bottom). As we see in section 4 and 
5, these structures are not specific to SIS models. By looking into these small-scale "fine 
structures" , one can measure the density profile and the mass of the lens perturber. 



3. SIS lens + background shear and convergence 

In this section, we show that in the large source size limit, the order of the (de)magnification 
perturbation of a circular top-hat source with radius L (in units of an Einstein radius) owing 
to an SIS subhalo is typically 0{L~ X ) with respect to that for an unperturbed lens. We 
model lensing by a subhalo locally as a background constant shear and convergence (Finch 
et al. 2002; Keeton 2003). If the mass scale of a substructure is sufficiently smaller than 
that of the macrolens, this model provides a good approximation in representing the local 
property of lensing along the line of sight to a substructure. For simplicity, we study a 
circular top-hat source centered at the SIS lens center. 



3.1. Systematic distortion 



First, we study the systematic distortion of a circular source. Let us consider an SIS 
lens with a one-dimensional velocity dispersion a and an external convergence k and shear 
7 > 0. In the coordinates aligned to the shear, the lens equation is 

v = (i-r)0 a " IJ ' 



where 



and 



r = 



K + 7 
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k — 7 



Br, 



6 E = 4tt 



a 2 D 



LS 



(13) 



(14) 



(15) 



c* D s 

where Ds is the angular diameter distance to the SIS lens, D^s is the distance between the 
SIS and the source, and c is the light velocity. In terms of normalized coordinates, x = 6 x /8e 
and y = 6 y /0E, the lens equation can be reduced to 



fi-rw- 



X 

Ixl 



(16) 
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To characterize the effect of background shear and convergence, we examine how a set of 
orthogonal vectors (y±, 0), (0, y 2 ) in the source plane is mapped to another set of orthogonal 
vectors (xi,0), (0,x 2 ) in the lens plane that is aligned to the shear (i.e., astrometric shifts). 
From the lens equation (16), we have 



On the one hand, because of equation (17), the presence of an SIS perturber leads to a shift 
of a horizontal vector of the unperturbed image (x± > 1,0) to (xi + 1/(1 — k — 7), 0), and 
a shift of a vertical vector of the unperturbed image (0,x 2 > 1) to (0,x 2 + 1/(1 — k + 7)). 
On the other hand, as we have seen in section 2, for k — 7 = 0, a circular top-hat source is 
isotropically magnified to a circular top-hat image with radius 1 + L. From these properties, 
we can expect that the boundary of the perturbed macrolensed image takes the form of an 
elliptic with a semi-major axis a+ 1/(1 — k — 7) and a semi-minor axis b + 1/(1 — k — 7). 
The validity of this approximation is discussed in section 3.2. 

Based on this approximation, distortion patterns are classified into three regimes by the 
sign of the eigenvalues of 1 — T, £1 = 1 — (« + 7) and £ 2 = 1 — (k, — 7):(1) £1 > 0,^2 > 
positive parity; (2) £1 < 0,^2 > negative parity; and (3) £1 < 0,^2 < doubly negative 
parity 1 . As shown in Figure 5, in the positive parity case, both the semi-major axis and 
the semi-minor axis of an image extend by Aa = £f 1 and Ab = l;^ 1 in comparison with 
the unperturbed lensed image (background shear and convergence only) while in the doubly 
negative parity case, both axes shrink, i.e., Aa < and Ab < 0. In the negative parity case, 
the semi- major axis shrinks by \Aa < 0| = l^f 1 ) whereas the semi- minor axis extends by 
Ab = £ 2 ~ 1 (Figure 6). Note, however, that a deviation from an elliptic is no longer negligible 
in the negative parity case, which is discussed in section 3.2. 



In this subsection, we examine systematic (de) magnification for a lens that consists of 
an SIS plus a constant shear and a convergence that depends on the sign of the eigenvalues 
of 1 — r. For simplicity, first, we consider a circularly symmetric top-hat source with radius 



1 These parities correspond to a minimum (positive), a saddle (negative), and a maximum (doubly nega- 
tive) point in the arrival time surface. 





3.2. Systematic (de) magnification 
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L > 1 centered at the position of an SIS. Next, we compare our analytically calculated 
(de)magnification values with numerically calculated ones in cases where the source is placed 
at an off-center position. 

First, we study the magnification effects in the positive parity case and the doubly 
negative parity case. As shown in Appendices B and C, in these cases, the shape of a distorted 
image perturbed by an SIS is well approximated by an ellipse, provided that L ^> I and 
1 1 — Avj —7 ^> 0. Therefore, a (de) magnification ratio defined as the ratio of the magnification 
for a perturbed lens (perturber + macrolens) to the magnification for an unperturbed caustics 
lens (macrolens only) x(L) = /V/ib g can be simply written as 

x(i)Kl + « + _L (18) 

Thus, the order of the leading term in the (de) magnification perturbation defined as x ~~ 1 
is ~ 0(L~ X \ A perturbed image is always magnified (demagnified) in comparison with the 
unperturbed image in the positive (doubly negative) parity case. As shown in Figure 7 (left), 
the analytic formula (18) for a circular source centered at an SIS is accurate to within a few 
percent for L 1. Even if a circular source is not centered at the SIS perturber, formula 
(18) still gives good accuracy for L ^> 1 as shown in Figure 8. 

Next we consider the negative parity case. If the shape of an image perturbed by an SIS 
is well approximated by an ellipse, the (de) magnification ratio should be x{L) = 1 — 1/L 2 . 
However, in the negative parity case, distortion from an elliptic shape is no longer negligible 
at 0(L~ r ). To quantify departure from an ellipse, we introduce two parameters e{9) and 
rj(9) in the lens equation in polar coordinates that satisfy 

Lcos(9 + rj) = a(l + e)(l — k — 7) cos# — acos9/R[a, b, 9} 

Lsm(9 + r]) = 6(1 + e)(l - k + 7) sin 9 - fcsin 9/R[a, b, 9], (19) 

where 

a = (L + sgn(£i))£r\ b = (L + sgn(£ 2 ))^\ (20) 

and 

R[a, b, 9] = (a 2 cos 2 9 + b 2 sin 2 9) 1/2 . (21) 

In terms of e, the ratio of the perturbed magnification factor to the unperturbed mag- 
nification factor can be approximately written as (see Appendix B) 

4 r /2 

X~~ d#((l + e)cos# + e'sin#)(l + e)cos#, (22) 
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where e' = de/dO. As shown in Appendix B, the order of e is Thus the order of 

the (de)magnification perturbation x{L) — 1 is again 0(L~ l ) in the large source size limit 
L >> 1. Note that the analytic formula (22) can be used for any type of parity. 

As shown in Figure 7, for L > 10, the analytic formulae (18) and (22) for a source 
centered at an SIS give an accuracy to within a few percent with respect to the numerically 
calculated values. Because we expanded x i n L" 1 , the expected error at the order 0(L~ V ) 
is Ax ~ C(L" 2 )/(1 + OiL' 1 )) ~ (L(L + For instance, we expect an error A% ~ 0.2 

for L = 2 but Ax ~ 0.01 for L = 10. If the distance between the boundaries of the source 
and caustics is very small, the magnification ratio x is sensitive to the position of the lens. 
However, in the large source size limit L 1, x asymptotically converges to values given by 
formulae (18) and (22). In this shown in Figure 8, (de) magnification is less sensitive 

to the position of the source provided that the caustics are totally occulted by the source. 
Our analytic formulae are very useful in representing the effects of background parameters 
k and 7 on systematic (de) magnification in such a case. 

As we have seen, for typical values of convergence k and shear 7, the (de) magnification 
perturbation 1 — x{L) can be as large as 1 — x{L) ~ 0.4 at L ~ 5 in the positive or doubly 
negative parity case. In other words, even if the source size is larger than the Einstein ring 
or caustics of an SIS perturber, the (de) magnificat ion effects owing to the perturber are still 
prominent. In the negative parity case, the (de) magnification perturbation is smaller than 
for other parities 1 — x{L) ~ 0.1 at L ~ 5 for typical values of k and 7 but it still cannot be 
negligible at the order C(L _1 ). 



3.3. Mean (de) magnification for a point source 

Our results for the systematic (de)magnification effects for a circular top-hat source 
centered at the lens center can be applied to the study of the mean (de)magnification for a 
point source within a radius L from the lens center. 

For a point source at y = y s in the source plane in which a lens perturber is put at the 
center, the magnification factor corresponding to an image x = x(y s ) is written as 



H = 



<9x 



dy 



(23) 



y=y s 



Then the magnification factor averaged over a disk with radius L centered at the center in 
the source plane is 



Mi 



4vrL 2 y |y| < L 



(24) 
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This is equivalent to the magnification factor for a top-hat circular source with radius L 
centered at the lens center. 

The obtained mean magnification factor (fi) L is related to the lensing cross section a (5) 
which is defined as the cross section for a magnification perturbation stronger than 8 (Keeton 
2003). Conversely, we can define the mean magnification perturbation for a cross section 
a = 4ttL 2 centered at the perturber as (6)l = 1 — (/-<■) L/Vbg- The mean magnification is 
useful for lensing systems in which the macro-lensing parameters are precisely determined. 
Comparing an observed value of the perturbation 5 with a computed mean perturbation 
and its variance, one can constrain the size of the cross section in units of the area of an 
Einstein disk of a perturber, provided that the source size is negligible in comparison with 
the Einstein ring. 

Thus, our results are useful not only for estimating the magnification for extended 
sources larger than the Einstein ring but also for estimating the mean magnification per- 
turbation for a point source in substructure lensing (see also Metcalf & Madau 2001; Chiba 
2002; Dalai & Kochanek 2002). 



4. Measurement of velocity dispersion and distance 

In this section, we study how one can break the degeneracy between the lensing effects 
of the velocity dispersion (or total mass) of a lens perturber and those of its distance along 
the line of sight by measuring distortion in the extended images. 

Suppose that a halo at redshift z — z^ deflects light from a source at redshift z s by an 
angle ach and a clump at redshift z — z c < Zh deflects the light by an angle a c (see Figure 
9). The normalized lens equations for an unperturbed system and a perturbed system are 

y = x - a h (xo), (25) 

and 

y = x - a c (x-x c ) - a h (x'), (26) 

respectively (Schneider et al. 1992). Here, y is the position of a source, x is the position of 
an image in the microlens plane, x c is the position of a clump, and x' is the position of an 
intersection between the light ray and the macrolens plane. The normalized angle is defined 
as ct x = D(z x , z s )/ D s 6e, where D(z x , z s ) is the angular diameter distance between an object 
at redshift z x and a source at redshift z s , D s is the distance to the source, and 0e is the 
Einstein radius of the macrolens. 



- 12 - 



From the lens equation for the microlens system, we have 

x' = x-/3a c (x-x c ), (27) 

where the factor 

= (28) 

D cs D h 1 ; 

encodes a distance ratio written in terms of the distance between the clump and halo D c h, 
the distance between the clump and source D cs , and the distance to the halo Dh- In what 
follows, we assume that the Einstein radius of the clump is much smaller than that of the 
halo, i.e., |x' — x | <C 1. Then the deflection angle is approximately given by (Keeton 2003) 



tth(x') ~ a/,(xo) + T(x - /3o! c (x - x c ) - x ), T = 



dct h 



<9x 



(29) 



X() 



where T represents the convergence and shear for the unperturbed macrolens. From equa- 
tions (27), (28) and (29), we obtain an effective lens equation for a perturbed system, 

y = (l-r efr )x-a c (x), (30) 

where 

x = x - x c , y = -(1 - r)(l - PT)-\ Xc - x ), (31) 

and 

r cff = i-(i-r)(i-/?r)- 1 . (32) 

Here y is the effective source position that satisfies y = (1 — /3r) _1 y if the convergence k and 
the shear 7 are constant. As a function of the convergence k and shear 7 in the unperturbed 
macrolens and the distance ratio f3, the effective convergence and shear can be expressed as 

(1 -/?)(/, -/3(/. 2 - 7 2 )) , (l-/?)7 r „x 

cS (1 - /3«) 2 - /3 2 7 2 ' 7cff (l-/3«) 2 -/3 2 7 2 ' 

In what follows, we assume that T is a constant matrix that does not depend on the image 
position x , i.e., = T, and we also assume that an SIS perturber is placed at the center 
of the image plane x c = (0,0) , i.e. ck c (x) = 8 Ep x/|x|, where 9 Ep is the Einstein radius of 
the SIS perturber. 

Now consider an astrometric shift of an image Ax = x — x that is defined as the dif- 
ference between a perturbed macrolensed image position x and an unperturbed macrolensed 
image position x that satisfies y = (1 — r)x . From equation (31), one can see that an effec- 
tive image position x that satisfies y = (1 — r cff )x is equal to the unperturbed macrolensed 
image position x . Then the astrometric shift of an image with respect to the unperturbed 
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image can be written as Ax = x — x . Taking coordinates x = (x±, x 2 ) aligned to the shear, 
astrometric shifts of images on the horizontal axis x\ are Axy = #e p /(1 — k c h — 7 e g) and those 
of images on the vertical axis x 2 are Axj_ = #e p /(1 — « e ff + Teff)- Plugging the relations in 
equation (33) into these shifts, one can see that the distance ratio j3, and the Einstein radius 
6e p of the perturber can be determined from observed values of astrometric shifts Axy and 
Axl provided that parameters k and 7 are determined from the best-fit macrolens model. 
Furthermore, if the macrolens parameters such as D s , D h , and D hs are known, then one can 
measure the distance between the SIS perturber and the macrolens from the distance ratio 
P and the velocity dispersion of the SIS. 

The accuracy in the measured distance ratio (3 for a given resolution can be estimated 
from an astrometric shift Ax = x — x of a perturbed image with respect to the unper- 
turbed image. As we have seen, in the coordinates aligned to the background shear, the 
f3 dependence of the horizontal shift for unperturbed images on the horizontal axis xi is 
Ax\ = 9e p (1 — P(k + 7))/(l — k — 7)- As (3 increases from to 1, A^i decreases by 
A(Axi) = 0E P (/c+7)/(l — k — 7). Therefore, for macrolenses with (/c+7)/(l — k — 7) = 0(1), 
the necessary angular resolution is A9 ~ 9e p in order to determine (3 with a relative accu- 
racy A/3//3 < 1. In other words, (3 can be measured from a horizontal shift A27 with good 
accuracy if the angular resolution is comparable to the Einstein radius of the SIS perturber. 
Similarly, the (3 dependence of the vertical shift for unperturbed images on the vertical axis 
x 2 is Ax 2 = 9e p (1 — (3(k — 7))/(l — /t + 7). Therefore, for macrolenses with k ~ 7, the (3 de- 
pendence of vertical shifts Ax 2 can be observed with an angular resolution A9 ~ (« — 7)#e p - 
In other words, one needs to resolve an image that is significantly smaller than the Einstein 
radius of an SIS perturber if k ~ 7. However, as we have seen in section 2, the Einstein 
radius 9e p can be directly measured by the size of the dipole structure where the surface 
brightness gradient at the position of the SIS perturber is non vanishing (see also section 5). 
Therefore, we conclude that (3 can be measured from the observed horizontal astrometric 
shift Axy with respect to an unperturbed image if observation with an angular resolution 
A9 ~ 9e p is achieved. 

In practice, however, we should pay attention to various factors that complicate the 
lens system. First, the observed distortion can be caused by distortion of the source itself. 
However, this is not a serious problem for multiply imaged lensing systems. By comparing 
a lensed image with other images, we can easily distinguish whether the distortion pattern 
is associated with the source properties or not. Second, the SIS perturber might not be 
spherically symmetric. Third, our assumption of an SIS density profile for perturbers may 
not be correct. Fourth, perturbers may not be single. Finally, the center of the perturber 
may not lie on the center of the source. The second and the third problems can be solved 
by measuring the distortion inside the Einstein radius of the perturber, which is discussed 
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in section 5. To solve the fourth and the last problems, we need a source with spatially 
varying surface brightness. Substructures in a QSO jet are good examples. If the position 
of the perturber is known, we can draw a set of coordinates for which the perturber resides 
at the center. Then we can assign position vectors to each substructure in the extended 
source, which will restrict the lens parameters that can characterize the microlens such as (3 
and 9e p - In other words, a sufficient number of substructures of an image can restrict the 
substructure lensing parameters. 



5. Measurement of mass density profile 

From a reconstruction of a mapping between a macrolensed source image and a mi- 
crolensed image, we can extract information of the mass density profile of the perturbers. 
To do so, we need to resolve distorted images within the Einstein radius of a perturber. Let 
us first consider a simple lens model that consists of a source S with a circular symmetry 
whose surface brightness obeys a Gaussian distribution with a standard deviation a, 



fg(y>°) = -/^ ex P 



y 2 



2a 2 



(34) 



and a lens (i.e. a perturber) P with a circular symmetry. For simplicity, we consider a 
polynomially suppressed mass density profile for the lens P. 

First, we consider a system in which a lens is put in the direction of the line of sight to 
the center of the source. Then we have the lens equation 

y = x ~R^' (35) 

where a corresponds to the power of the polynomially suppressed mass density profile of 
P. For instance, a = 2 corresponds to a point mass, and a — 1 corresponds to an SIS. 
From equation (34) and (35), for a > 0, we obtain the surface brightness profile f(x) in the 
neighborhood of a lens x = \x\ <C 1 as, 

In/(x)«-^p-lnv^F(7. (36) 

In Figure 10, one can clearly see the difference in the surface brightness of lensed images that 
depends on the mass density profile of the lens. For a lens with 1 < a, the surface brightness 
vanishes toward the center of a lens as ln/(x) ~ — x 2 ~ 2a . This is because the outer dimmed 
region of the source is mapped into the neighborhood of the center of the lens. Consequently, 
one would observe a "black hole" at the position of a point mass lens if a — 2. In contrast, 



- 15 - 



the surface brightness is finite in the neighborhood of the center of an SIS lens. In this case, 
a region R : \y\ < 1 is mapped into a region R! : \x\ < 1 on the opposite side (with opposite 
sign) within the Einstein radius but any regions outside the Einstein radius are not mapped 
into a region inside the Einstein radius. Hence, one would see a "dark hole" instead of a 
"black hole" for an SIS lens. If the mass density profile of a lens is shallower than that of an 
SIS lens, i.e. a < 1, then one would see a bright spot surrounded by a dark hole, because 
the mass concentration is so weak that the bending angle for a light ray that departed from 
the brightest region of the source is very small. 

In real observations, the detectability of these features strongly depends on the resolution 
of the map and the linear scale a in which the surface brightness varies. As a gets larger, 
variability in the surface brightness of the lensed image within an Einstein radius gets smaller 
as oc a^ 1 (see Figures 10 and 11). To determine a or equivalently, the Einstein radius of 
a perturber and a, one needs to measure a mapping between the source plane and the 
lens plane. Suppose that we observe a "dark hole" at a peak in the surface brightness of 
the source. If there are peaks or dips at yi > 1 in the source, then mapped images of 
these structures should be observed inside dark holes at Xi < 1. From more than two sample 
points of these structures, one can determine the model parameters a and a. Even if circular 
symmetry in the projected mass distribution of the perturber is broken, a sufficient number 
of sample points can determine the lens parameters such as ellipticity and external shear that 
describe distortion from a circular symmetry. For instance, in a QSO-galaxy lens system, 
observation of astrometric shifts of subjects in a QSO jet with respect to a macrolensed 
image may determine the model parameters for perturber lenses. 

Next we consider lens systems in which the lens is not placed at the peak of the source. 
As shown in Figure 11, we would observe a pair of dark and bright spots at the position of 
the lens where the radial derivative in the surface brightness does not vanish. For 1 < a < 2, 
the surface brightness vanishes at the position of the lens because the outer dimmed region 
is mapped inside the Einstein radius. For a — 1, the surface brightness is finite at the 
position of the lens but the radial derivative of the surface brightness at the position of the 
lens diverges even if the surface brightness of the source is smooth. As in the former case, 
from more than two sample points, one can determine the model parameters o and a. 

Our assumption of a Gaussian form of the surface brightness of a source may be too 
idealized. In order to estimate the effects of deviation from the Gaussian distribution, we 
consider small-scale fluctuations in the background source. For simplicity, we added Gaussian 
fluctuations with a standard deviation that is 1/15 of the background Gaussian peak value 
and with a vanishing mean to the smooth background Gaussian source. We calculated lensed 
images for two examples with the same smoothed source parameters as in the previous cases, 
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namely, a — 2 and y = — 1 (large surface brightness gradient), and a — 10 and y = — 5 (small 
surface brightness gradient). The minimum wavelength of the Fourier mode of an additional 
fluctuation is assumed to be A m i n = 1/4. As shown in Figure 12, a discontinuous change 
in the surface brightness for a point mass is still apparent for both cases, whereas for an 
SIS, such a discontinuous change is less apparent for the case with a small-scale brightness 
gradient (right). This is because for an SIS, an observable astrometric shift is limited to 
images within the Einstein ring of the perturber. Therefore, in order to measure the density 
profile of an SIS, we need a source whose (smoothed) surface brightness gradient is sufficiently 
large. 



6. Summary 

In this paper, we have explored the extended source size effects in substructure lensing. 
First, we have shown that for a simple SIS lens model with a circular top-hat or Gaussian 
surface brightness profile, the magnification effect is prominent even if the source size is 
larger than the size of an Einstein ring. Second, we have analyzed SIS lens systems with 
background shear 7 and convergence k relevant to substructure lensing. We have derived 
analytic (de) magnification asymptotic formulae in the large source size limit for three types 
of parity of an image, namely, positive, negative, and doubly negative parities. In the positive 
(doubly negative) parity case, the source is always magnified (demagnified) in comparison 
with the unperturbed macrolens. In the negative parity case, magnification depends on 
the sign of 1 — k. For 1 — k > 0, the source tends to be slightly magnified, whereas it 
tends to be slightly demagnified for 1 — k < 0. Again, we find that in the large source 
size limit, the order of the (de) magnification perturbation 1 — x(L) with respect to the 
unperturbed one is typically C(L _1 ), where L is the source radius in units of an Einstein 
radius. We have also shown that these results are relevant to the mean (de)magnification 
for a point source at a distance L from the lens perturber center. Third, we have shown 
that one can break the degeneracy between the lensing effects of mass and distance of a 
substructure based on the distortion pattern in the perturbed image provided that the lensing 
parameters of macrolensing are determined from the position and flux of multiple images. 
Finally, we have shown that the density profile of a substructure is directly measured by 
reconstructing the mapping between a perturbed lens system (macrolensing + microlensing) 
and an unperturbed one (macrolensing only), which can be achieved by resolving an image 
within the Einstein radius of a substructure. 

Although it seems that our results may not be relevant to current observations of sub- 
structure lenses, they will surely become important tools in future observations with im- 
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proved sensitivity and angular resolution. For instance, if the sensitivity in the measured 
flux is dramatically improved, then we may be able to find an asymptotic regime in which 
the magnification perturbation is inversely proportional to the size of the source with a com- 
mon center of light. Future observations of QSO jets in the radio band or of cold dusts in 
starburst galaxies in submillimeter band may reveal such systems if observed with different 
frequencies. If the proper source size is known, then we can determine the angular size of 
the Einstein radius of a lens perturber. Similarly, if the angular resolution were dramatically 
improved, then we would be able to get a fruitful amount of information from a substructure 
within a source. For instance, a QSO jet usually consists of two jet components plus a center 
nucleus. Each jet sometimes has some small subjects or small-scale structures. Then we 
could determine the spatial variation within the source in the global linear maps between 
macro lensed images by comparing their substructures. Mapping cold dusts around a QSO 
in a QSO-galaxy lens at the submillimeter band with an angular resolution of ~ 10 mas may 
directly determine the distance to a ~ 10 s — 10 9 M Q subhalo near the center of the lens galaxy 
(Inoue & Chiba 2005). Breaking the mass-distance degeneracy and a precise measurement 
of mass density profile of a perturber might both be achieved by such an observation. 

In reality, lensing systems are usually much more complex and modeled with a large 
number of parameters (often too many!). For instance, a macrolens can be perturbed by 
a number of subhalos with different projected masses, ellipticities, and distances from the 
source. The corresponding best-fit parameters can be determined by minimizing a certain 
measure between an observed image and a reconstructed image based on a certain lensing 
model. 

Many important issues that we have omitted are (1) irregularity in the source shape, 
(2) irregularity in the mass density profile of each perturber, (3) time varying source effects, 
(4) the effect of multiple perturbers, and (5) lensing parameter error estimates. In real 
observations, we always have to tackle all these problems. We will soon address these issues 
in our future work. 

We would like to thank S. Kameno, T. Minezaki, N. Kashikawa, and H. Kataza for 
useful discussions on observational aspects of substructure lensing. We would also like to 
thank the anonymous referee for useful comments. This work has been supported in part 
by a Grant-in-Aid for Scientific Research (15540241) of the Ministry of Education, Culture, 
Sports, Science and Technology in Japan. 



A. Analytical treatment of overlapping region 



We consider a circularly symmetric top-hat source disk D with a radius L located at 
a distance Co from the center of an SIS lens. The magnification factors for an image with 
positive parity and negative parity can be written in terms of the radius L and Co as in 
equation (5). 

For Co ~ L > 1, there is no overlapping region between a source disk D and a unit 
Einstein ring E, so the total magnification is // = //+ + //_ = For ( — L < —1, an 
Einstein ring E is totally included in a source disk D, and so we have // = pL + + 7r/z_(L = 
1, Co — 0)/tcL 2 = + 1/L 2 . For Co < 1, if a source disk D is totally included in an Einstein 
ring E, i.e., L < |1 — £o|, then fj, — /i + + //_. These results are summarized in table 1. 

For other cases in which there are some overlaps between a source disk D and an 
Einstein ring E, namely, |Co — 1| < L < ( + 1, it is difficult to obtain a simple explicit 
formula. However, for the cases L > 1 and Co 1 overlapping some part of an Einstein 
ring E, a simple approximated formula can be obtained as follows. In coordinates (L, fJ>(L)), 
we may connect two points A, (Co — l,//+(Co — 1))) an d B, (Co + l ! / i +(Co + 1)), simply by 
a segment, because the area of the overlapping region E U D is a monotonically increasing 
function of L. However, the derivatives at A and B are not continuous. Instead, we can 
connect A and B by a smooth function whose derivatives vanish at A and B, because the 
derivatives d/i/dL at L = Co — 1 and at L = Co + 1 are very small for Co ^ 1- A simple choice 
of such a function is 

x,,(*Tf\ Erf( S (L-Co))-Erf(- S (L-Co)) + 2 

<W> L > Co) = ^2 ' ( A1 ) 

where s > controls smoothness of 5ji around L = Co- Then we can approximate the 
magnification factor /j, as 

MA Co) Co) + Co), L>l,Co>l- (A2) 

Table 1. Three Different Cases for 
Magnification of a Circular Top-hat Source 3, 



Case no overlapping D C E E C D 



a Here D and E are a source disk and a unit 
Einstein ring, respectively. The definition of /i + 
and /i_ is given in equation (5). 
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We find that s ~ 3 leads to an accuracy within a few percent for Co > 1 (see Figure 1). Note 
that our approximated formula (A2) fails for Co < 1 and L < 1 + Co because the contribution 
from ji_ cannot be negligible in that case. If Co ^ 1> then /i can be approximated as 

H(L, Co) ~ H+(L, 0) + n_(L, 0)H(1 - L), (A3) 

where H is a Heaviside step function. 



B. Deformation parameters in the negative parity case 

Assuming 77 <C 1, equation (19) can be solved as 

e (0) = (6(i_ K + 7 ))- 1 jL-6((l-/t + 7 )-ir 1 ) - ^cos 2 9 (b(l - K + -f)LR 

+ a(26 7 — (1 — k — i)LR)) ) /R(a(-1 + k + 7) cos 2 + 6(-l + k - 7) sin 2 0) tel) 



and 



W = 1 2flx« -2/) • ( B2 ) 

a£i cos z y + 0^2 sin 6* 



In the negative parity case, the approximated solution of e ~ 0(L x ) for L ^> 1 is 
e(0) « L~ 



■?n sin 2 cos 2 9 

-1 + 2 cos 2 + — — + 



(1 - « + 7 )F(«, 7, 0) (1 - « - 7 )F(«, 7, 0) 



(B3) 



where 



cos 2 # sin 2 # 
^■7,^) = ^ ;i _ K _ 7)2 + (1 _ K + 7)2 . (B4) 

Now, we evaluate a mean value of e averaged over 9 in the negative parity case. For 1 — k > 0, 
we have |1 — /t + 7 |~ 1 < |1 — k — 7 | _1 . Therefore, the contribution of the integrand cos 2 9/(1 — 
k — 7)F is dominant over sin 2 9/(1 — k + 7 )F. Note that the amplitudes of both integrands 
are 1. Because (—1 + 2 cos 2 9} e = 0, we have (e) g < 0. Thus, the area of the perturbed image 
tends to become smaller than that of an elliptic disk D with a semi-major axis a and a 
semi-minor axis b at 0(L~ V ). For 1 — k < 0, we have |1 — k + 7 |~ 1 > |1 — k — 7 |~ 1 . Therefore, 
the contribution of the integrand sin 2 0/(1 — k + j)F is dominant over cos 2 0/(1 — k — 7) F, 
leading to (e)$ > 0. In this case, the area of the perturbed image tends to become larger than 
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that of E (see Figure 6). In summary, we expect a magnification ratio x > 1 for 1 — k < 
but demagnification, x < 1 for 1 - k > 0. 

As for ?7, equation (B2) yields 

77(0) « L _1 sin20( 1 + xo 7 0Xr ^ -V (B5) 

Because 7/((l — k) 2 — 7 2 )-F < 1, we have < L^ 1 . Thus, an assumption of rj <C 1 is correct 
for L > 1. 



C. Deformation parameters in the positive and doubly negative parity cases 

To evaluate the accuracy of equation (18), we calculate the deformation parameters e 
and r\. In the positive parity case, the semi-major axis and semi-minor axis of an ellipse D 
that represents an approximated shape of the perturbed lensed image are a = (L + l)£j~ 1 
and b = (L + l)^ 1 . Plugging these values into equation (B3), we obtain 

,(ft\ t-i( l-/t + 7 cos2fl \ 

The order of e is evaluated in terms of the median (1 — k)/((1 — k) 2 — 7 2 ) > of F as 

eW ~L-.f2f^Y ( C2 ) 

Therefore, for 1 — k — 7 ^> 0, the effect of deformation from an ellipse D is negligible at 
O(L^). 

In the doubly negative parity case, we have a = (L — 1 and b = (L — l)^ 1 . In a 
similar manner, one can show that 

Therefore, for — 1 + k — 7 3> 0, the effect of deformation from an ellipse D is again negligible 
at OiL- 1 ). 
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Fig. 1. — Magnification factor \x (solid curves) for a top-hat circular source lensed by an 
SIS without background as a function of a source radius L with distance from lens center 
Co = 1 (top curves), Co = 2 (middle curves) and Co = 3 (bottom curves) in comparison with 
analytically calculated values based on eq. (6) with smoothing parameter s = 3 (dashed 
curves). 
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Fig. 2. — Lensed images of a top-hat circular source centered at ( — (3, 0) with a radius 
L — 1 (top left), L = 2 (top middle), L = 2.5 (top right), L = 3 (bottom left), L = 3.5 
(bottom middle), and L = 4 (bottom right) in the source plane. An SIS lens is centered at 
(0, 0) in the lens plane. The Einstein radius is set to unity. 
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Fig. 3. — Magnification factor ji {solid curves) for a Gaussian circular source lensed by an 
SIS without background as a function of a standard deviation L with distance from the lens 
center Co = 1 {top curves), Co = 2 {middle curves), and Co = 3 {bottom curves) in comparison 
with those for an image with positive parity only (see eq. (10); dashed curves) and those for 
top-hat sources with a scaling L — > a/8/7tL (here L denotes a radius of a circular source) 
for the same Co values {dotted curves). 
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Fig. 4. — Top: Lensed images of a Gaussian circular source centered at £ = (3, 0) with a 
standard deviation T = 1 (Ze/t), L = 2 (middle), and T = 3 (right). Bottom: Lensed images 
with the same parameters except for the position of the source at £ — (0,0). An SIS lens 
is centered at (0,0). Each interval between neighboring contours is 1/10 of the maximum 
value in the surface brightness. The Einstein radius is set to unity. 
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Fig. 5. — Systematic distortion by an SIS perturber at the center of a top-hat circular 
source with radius L ^> 1 in an external shear 7 and convergence k. The dotted and solid 
closed curves show the boundaries of an unperturbed and a perturbed macrolensed image, 
respectively. 
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Fig. 6. — Systematic distortion by an SIS perturber at the center of a top-hat circular 
source in an external shear 7 and convergence k. The dotted and solid closed curves show 
the boundaries of an unperturbed and a perturbed macrolensed image, respectively. The 
dash dotted closed curves are boundaries of an ellipse with a semi-major axis a = (L — l^f 1 
and a semi-minor axis b = (L + l)^ 1 . The total area inside the boundary gets smaller 
(larger) than the unperturbed macrolensed image for 1 — k > (for 1 — k < 0) provided that 
L > 1. 
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Fig. 7. — Comparison between numerically calculated magnification ratios \ {solid curves) 
and analytically calculated ones (dashed curves). For the maximum (positive parity) and 
minimum (doubly negative parity), eq. (18) is used, and for saddles (negative parity), eqs. 
(Bl) and (22) are used. Parameters are (/c, 7) =(0.3,0.3) (top left), (2.5,0.2) (bottom left), 
(0.6,0.6) (top right), and (1.8,2.0) (bottom right). The center of an image is placed at y =(0,0) 
in the source plane. The Einstein radius of an SIS lens is set to unity. 
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Fig. 8. — Magnification ratios x as a function of radius L of a circular top-hat source for 
(«, 7 ) =(0.3,0.3) (top left), (2.5,0.2) (bottom left), (0.7,0.7) (top right), and (1.8,2.0) (bottom 
right) for different image positions centered at y =(2,0) (solid curves), (1,1) (dashed curves), 
(0,2) (dash dotted curves) in the source plane. The Einstein radius of an SIS lens is set to 
unity. An increase at L < 2 in the minimum and saddle cases is owing to a crossing with a 
caustic. 




Fig. 9. — Light from a source deflected by a halo at Zh and a clump at z c < Zh- 




Fig. fO. — Left: Surface brightness of lensed images for a = 2 (thin solid curve), 1.5 (dashed 
curve), 1 (dotted curve), and 0.8 (dash dotted curve) for a circularly symmetric Gaussian 
source with a = 1 (thick solid curve). The horizontal axis denotes a coordinate x in the 
image plane. The Einstein radius of a perturber is set to unity. A lens is put at the peak 
of the Gaussian distribution. Right: Same parameters as in the left panel, expect for the 
standard deviation a = 4. 
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Fig. 11. — Same parameters as in Figure 10 expect for the standard deviation a and the 
position of the lens. Left, a = 2 and a source is put at y = — 1 in the source plane. Right 
a — 10 and a source is put at y = — 5 in the source plane. 
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Fig. 12. — Effects of small-scale fluctuations in the surface brightness of the source. The 
parameters are the same as in Fig. 11, except for additional Gaussian fluctuations with a 
standard deviation equal to 1/15 of the Gaussian peak of the smooth background and a 
mean equal to 0. The minimum wavelength of the fluctuation is assumed to be A m i n = 1/4. 
Thick solid curves denote the surface brightness / of the source, and thin solid curves and 
dotted curves correspond to the surface brightness of an image for a = 2 (point mass) and 
a = 1 (SIS), respectively. 



